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Abstract 

Scattering of electromagnetic (EM) waves by one and many small 
(ka <g; 1) impedance particles Dm of an arbitrary shape, embedded 
in a homogeneous medium, is studied. Analytic formula for the field, 
scattered by one particle, is derived. The scattered field is of the order 
0(0^""), where k e [0, 1) is a number. This field is much larger than in 
the Rayleigh-type scattering. An equation is derived for the effective 
EM field scattered by many small impedance particles distributed in a 
bounded domain. Novel physical effects in this domain are described 
and discussed. 
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1 Introduction 

In this paper we develop a theory of electromagnetic (EM) wave scattering by 
one and many small impedance particles (bodies) Dm, 1 < m < M = M{a), 
embedded in a homogeneous medium which is described by the constant 
permittivity eo > and permeability /ip > 0. The smallness of a particle 
means that ka <C 1, where a = 0.5diamDm is the characteristic dimension 
of a particle, k = ci;(eo/io)^^^ is the wave number in the medium exterior 
to the particles. Although scattering of EM waves by small bodies has 
a long history, going back to Rayleigh (1871), see [1], [16], the results of 
this paper are new and useful in applications because light scattering by 
colloidal particles in a solution, and light scattering by small dust particles 
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in the air are examples of the problems to which our theory is applicable. 
The Mie theory deals with scattering by a sphere, not necessarily small, 
and gives the solution to the scattering problem in terms of the series in 
spherical harmonics. If the sphere is small, ka <^ 1, then the first term in 
the Mie series yields the main part of the solution. Our theory is applicable 
only to small particles. They can be of arbitrary shapes. The solution to 
the scattering problem for one small particle of an arbitrary shape is given 
analytically. For many such particles the solution is reduced to solving a 
linear algebraic system. This system is not obtained by a discretization of 
some boundary integral equation, and it has a clear physical meaning. Its 
limiting form as o — > yields an integro-differential equation for the limiting 
effective field in the medium where the small particles are distributed. Wave 
scattering problem in the case of one body can be studied theoretically only 
in the limiting cases of a small body, ka <^ 1, or a large body, ka ^ 1. In 
the latter case the geometrical optics is applicable. This paper deals with 
the case ka <C 1. Rayleigh (1871) understood that the scattering by a small 
body is given mainly by the dipole radiation. For a small body of arbitrary 
shape this dipole radiation is determined by the polarization moment, which 
is defined by the polarizability tensor. For homogeneous bodies of arbitrary 
shapes analytical formulas, which allow one to calculate this tensor with any 
desired accuracy, were derived in [1^ . These bodies were assumed dielectric 
or conducting in [16j. Under the Rayleigh assumption the scattered field 
is proportional to a"^, to the volume of a small body. The physically novel 
feature of our theory is a conclusion that for small impedance particle with 
impedance C = h{x)a~'^, where h{x) is a given function, Reh{x) > 0, and 
K G [0, 1) is a given number, the scattered field is proportional to a^"**, that 
is, it is much larger than in the Rayleigh case, since « a^"" if a << 1 
and K > 0. 

In this paper wave scattering by small impedance particles is studied. 
Besides high intrinsic interest in this problem, the theory we develop al- 
lows one to get some physically interesting conclusions about the changes 
of the material properties of the medium in which many small particles are 
embedded. The results of this paper can be used to develop a method for 
creating materials with a desired refraction coefficient by embedding many 
small impedance particles into a given material. Such a theory has been 
developed by the author for scalar wave scattering, for example, acoustic 
wave scattering, in a series of papers [3]-|14j. The novel physical idea in this 
paper is to reduce solving the scattering problem to finding some constant 
pseudovector Q (see formula ([TH]) ). rather than a pseudovector function J 
(see formula (flOl) ) on the surface of the scatterer. The quantity Q is a pseu- 
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dovector which somewhat analogous to the total charge on the surface of the 
perfect conductor with the shape of Dm, while the function J is analogous to 
the surface charge density. We assume for simplicity that the impedance (" 
(see formula ([5])) is a constant given in (jl6p . A similar assumption appeared 
in paper [3], where scalar wave scattering theory was developed. The results 
of the theory in [3] was a recipe for creating materials with a desired refrac- 
tion coefficient in acoustics (see [l2], [13], [H]). The boundary impedance 
([16]) grows to infinity as a — )• 0, and one can pass to the limit in the equation 
for the effective (self-consistent) field in the medium, obtained by embed- 
ding many small impedance particles into a given medium. Such a theory is 
briefly summarized in paper [13], where the equation for the limiting field in 
the medium is derived. The aim of this paper is to develop a similar theory 
for EM wave scattering by many small impedance particles embedded in a 
given material. 

For EM wave scattering by one small body D of an arbitrary shape with 
an impedance boundary condition an analytic formula for the electromag- 
netic field in the region r := \x\ ^ a, is derived: 

E{x) = Eo{x) + [V—,Q], ry^a, g{x,y) := —, r = \x - y\, 

where £"0 is the incident field, which satisfies Maxwell's equations in the 
absence of the scatterer D, [A,B] = A x B is the cross product of two 
vectors, {Q,ej) = Q • ej is the dot product, {ej}^^^ is an orthonormal basis 
in M^, 

Q, := (Q, ej) = -^H,p(V x Eo{0))p, E := {I + a)r, (2) 

over the repeated index p summation is understood from 1 to 3, C is the 
boundary impedance, 1 5*1 is the surface area of the particle, the matrix Ejp 
is defined by the formula 



{I + a)(^6jp-\S\-^ J^Nj{s)Np{s)ds'^ := {I + a)Tjp, (3) 



where Nj{s) is the j—th component of the unit normal N{s) to the surface 
5" at a point s G S, pointing out of D, t := I — b, b is a matrix, bjp := 
\S\~^ Nj{s)Np{s)ds, k = ui{€ofio)^^'^ is the wave number, O € -D is the 
origin, / is the identity matrix, and a is a matrix. 



I + a:={I + l3) 



-1 



3 



where the matrix j3 is defined in ()29p (see below). 

By S"^ the unit sphere in M'^ will be denoted. The boundary S of the 
small body D is assumed smooth: it is sufficient to assume that in local 
coordinates the equation of S is given as X3 = (/>(xi,X2), where the function 
4> has first derivative satisfying a Holder condition. 

Briefly speaking, there are three basic novel results in this paper. The 
first result are formulas ([I])-® and equations (fSTT) and (f32]l . The second 
result is the reduction of the solution to many-body scattering problem to 
solving a linear algebraic system (LAS), see equations (jlTI) and (ll2]) . The 
third result is a derivation of the equation for the limiting effective (self- 
consistent) field in the medium in which many small impedance particles 
are embedded, see equation ()15|) . 

The scattering problem by one small body is formulated and studied in 
Section 2, the reduction of the solution to the many-body EM wave scatter- 
ing problem to the solution of a LAS is given in Section 3. Also in Section 3 
a derivation of the equation for the limiting effective field is given as a — >■ 
and M = M(a) —7- 00. In Section 4 the conclusions are formulated. 

In this paper we do not solve the boundary integral equation to which the 
scattering problem can be reduced in a standard approach, but find asymp- 
totically exact analytical expression for the pseudovector Q which defines the 
behavior of the scattered field at distances d» a. 

In fact, these distances d can be very small if a is sufficiently small, and 
d can be much less than the wavelength A = 

Therefore our theory is valid in the physical situations in which multiple 
scattering effects are dominant. 

2 EM wave scattering by one small impedance 
particle 

Let us use in this Section the following notations: D is a small body, D' := 
M"^ \ D, /c > is a wave number, /ca <C 1, a = 0.5diamD, k = A is the 
wavelength of the incident EM wave, fc^ = w^eo/^o, where u is the frequency 
and eo, Ho are constant permittivity and permeability of the medium. Our 
arguments remain valid if one assumes that the medium has a constant 
conductivity ctq > 0. In this case eo is replaced by eo + i^- Denote by 
S the boundary of D, by \S\ its surface area, by V the volume of D, by 
[E, H] = E X H the cross product of two vectors, and by {E, H) = E ■ H 
the dot product of two vectors, is the unit normal to S pointing out of 
D, ^ is the boundary impedance of the particle. 



4 



Let D be embedded in a homogeneous medium with constant parameters 
eo, /io- Electromagnetic (EM) wave scattering problem consists of finding 
vectors E and H satisfying Maxwell's equations: 



V y.E = ioj/doH, V X H = -iojeoE in D' ■.= R^\D 



(4) 



the impedance boundary condition: 



[N, [E,N]]=C[N,H] 



on S 



(5) 



and the radiation condition: 



E = Eo + VE, H = Hq + vh 



(6) 



where Eq,Hq are the incident fields satisfying equations ([4]) in all of M.^, 
ve = V and vh are the scattered fields. In the literature, for example in [I], 
the impedance boundary condition is written sometimes as -E* = A'^], 
where N is the unit normal on S pointed into D. Since our N is pointed 
out of D, our impedance boundary condition ([5]) is the same as in [1]. 

One often assumes that the incident wave is a plane wave, i.e., Eq = 
g^ika-x ^ £^ is a constant vector, a G 5^ is a unit vector, 5^ is the unit sphere 
in R^, a ■ 8 = Ve and vh satisfy the Sommerfeld radiation condition: 
—ikv) = o(l) as r := \x\ oo, and, consequently, [r^'jE^] = H + o{r^^) 
and \r^, H\ = —E + o(r~^) as r := |x| oo, where := x/r. 

It is assumed in this paper that the impedance C is a constant, Re ^ > 0. 

This assumption guarantees the uniqueness of the solution to Maxwell's 
equation satisfying the radiation condition. For completeness a proof of the 
uniqueness result is given in Lemma 1. The tangential component of E on 
S, £'*, is defined as: 



This definition differs from the one used often in the literature, namely, from 
the definition E^ = [N, E] . Our definition ([7]) corresponds to the geomet- 
rical meaning of the tangential component of E and, therefore, should be 
used. The impedance boundary condition is written in [Ij as E^ = C[H*', iVj], 
where C is the boundary impedance and Ni is the unit normal to S pointing 
into D. In our paper is the unit normal pointing out of D. Therefore, 
the impedance boundary condition in our paper is written as in equation 
([5]). If one uses definition ([7]), then this condition reduces to ([5]), because 
[[N, [H, N]],N] = [H, N]. The assumption ReC > is physically justified by 
the fact that this assumption guarantees the uniqueness of the solution to 
the boundary problem ([!])- ([6]). 



E^ = E - N{E, N) = [N, [E, N]]. 



(7) 
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Lemma 1. Problem (H])-® has at most one solution. 



Proof of Lemma [H In Lemma 1 it is not assumed that D is small. 
The proof is valid for an arbitrary finite domain D. To prove the lemma 
one assumes that Eq = Hq = and has to prove that then E = H = 0. 
Let the overbar stand for the complex conjugate. The radiation condition 
implies [-ff, x''] = E + o{\x\~^), where := x/r is the unit normal on 
the sphere Sr centered at the origin and of large radius R, and one has 
/' := fgjE,H] ■ x^ds = Jg^ \E\'^ds + o(l) asR^oo. 

From equations (jj]) one derives: 

/ (H -V X E - E -V xH)dx = i {iuj^iQ\H\^ - iujeo\E\^)dx, 
Jdr Jdr 

where Dr := D'nBji, and > is so large that D C Bji := {x : \x\ < R}. 
Recall that V ■ [E,H] = H- VxE — E-VxH. Applying the divergence 
theorem, using the radiation condition on the sphere Sr = OBr, and taking 
real part, one gets 

= -Re / [E,H] ■ Nds + Re / [E,H] ■ x^ds := I + l' . 
Js Jsr 

The radiation condition implies /' > as -R — > oo. The minus sign in front 
of the integral I comes from the assumption that on 5 is directed out of 
D. The impedance boundary condition and the assumption Re^ > implies 
/ > 0. One has I + I' = 0. One has 

- / [E,H] ■ Nds = j E- [N, H]ds = [ E- [N, [E,N]]/C, 
^ S J S J s 



and E ■ [N, [E,N]]/C = \E*\^C\Ct^- Therefore, / = ReCICrVs \Et?ds > 0, 
and /' > 0, the relation / + /' = implies / = 0. Therefore, if ReC > 
0, then Et = Q on S. Consequently, E = H = {) m. D. If ReC > 0, 
then lim/j^oo Jg \E\'^ds = 0, and since E satisfies the Sommerfeld radiation 
condition it follows that E = H = in D' . 

Lemma [1] is proved. □ 
Let us note that problem ([3])- ([6]) is equivalent to the problem 

W xV X E = k'^E in D', H = ^ ^ (8) 

[N, [E, N]] = [N, V xE]on S, (9) 
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together with the radiation condition ([6]). Thus, we have reduced the scat- 
tering problem to finding one vector E{x). If E(x) is found, then H = 
and the pair E and H solves Maxwell's equations, satisfies the impedance 
boundary condition and the radiation condition ([6]). 
Let us look for E of the form 

E = Eo + Vx g{x,t)J{t)dt, g{x,y) = ^ (10) 

Js 47r\x-y\ 

where £"0 is the incident field, which satisfies Maxwell's equations in the 
absence of the scatterer D, t is a point on the surface S, t € S, dt is an 
element of the area of S, and J{t) is an unknown pseudovector-function on 
S, which is tangential to S, i.e., N{t) ■ J{t) = 0, where N{t) is the unit 
normal to S at the point t € S. That J = J{t) is a pseudovector follows 
from the fact that £" is a vector and V x (gJ) is a vector only if J is a 
pseudovector, because g is a scalar. 

It is assumed that J is a smooth function on S, for example, J € C^{S). 

The right-hand side of (llOp solves equation ^ in D for any continuous 
J{t), because Eq solves ([8]) and 



VxVxVx / g{x,t)J{t)dt = graddivV X / g{x,t)J{t)dt 
Js Js 

-V^Vx / g{x,t)J{t)dt = k^V X / g{x,t)J{t)dt, x e D' . 
Js Js 



ill) 



Here we have used the known identity divcurlE = 0, valid for any smooth 
vector field E, and the known formula 

-\/'^g{x,y) = k'^g{x,y) + 5{x-y). (12) 

The integral Jg g{x,t)J{t)dt satisfies the radiation condition. Thus, formula 
pOj) solves problem ([8|), Q, ([6]), if J (t) is chosen so that boundary condition 
([9]) is satisfied. 

Let O € -D be a point inside D. The following known formula (see, for 
example, [2j) is useful: 



[iV,Vx / g{x,t)Jit)dt]^ = [ [Ns,[V^g{x,t)\^=s,J{t)]]dt± 
Js Js 



Jis) 



(13) 



where the it signs denote the limiting values of the left-hand side of (|13|) 
as X — 7> s from D, respectively, from D' . To derive an integral equation for 
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J = J{t), substitute E{x) from ([TO]) into impedance boundary condition ([9]), 
and get 

= /+[iV,[Vx / g(s,t)J{t)dt,N]]—^[N,VxVx [ g{s,t)J{t)dt], (14) 
Js *^/^o Js 

where 

/ := [N, [Eo{s),N]] - -^[AT, V X Eo]. (15) 
We assume that 

C = 4' (16) 

where Re h >0. and k G [0, 1) is a constant. 
Let us write pUj) as 

E(x) = Eo(x) + [V.5(2;, O), Q] + V X [ ig{x, t) - g{x, 0))J{t)dt, (17) 

Js 

where 

Q := / J{t)dt. (18) 

The central physical idea of the theory, developed in this paper, is: the 
third term in the right-hand side of (jl7p is negligible compared with the 
second term if fca ^ 1. Consequently, the scattering problem is solved if 
Q is found. A traditional approach requires finding an unknown function 
J{t), which is usually found numerically by the boundary integral equations 
(BIE) method. The reason for the third term in the right-hand side of (|17p 
to be negligible compared with the second one, is explained by the estimates, 
given below. In these estimates the smallness of the body is used essentially: 
even if one is in the far zone, i.e., ^ <C 1, one cannot conclude that estimate 
(pT]) (see below) holds unless one assumes that ka <^ 1. Thus, the second 
sum in pT|) cannot be neglected in the far zone if the condition fca <C 1 does 
not hold. 

We prove below that Q = 0{o?'~'^). To prove that the third term in the 
right-hand side of (jl7p is negligible compared with the second one, let us 
establish several estimates valid if a — )■ and d := \x — 0\ 3> a. Under these 
assumptions one has 

Ji := \[V.g[x,0),Q\\ < O fmaxjl^l^ ©(a^"''), (19) 



J2 := |V 



l^{g{x, t) - g{x, 0))a{t)dt\ < aO (max { ^}) 

(20) 
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and 

J2 



^ O (^max|^,A;a}) ^ 0, ^ = o{l), a ^ 0. (21) 
These estimates show that one may neglect the third term in (jl7p . and write 



E{x) = Eoix) + [V.gix, O), Q]. (22) 

The error of this formula tends to zero as a ^ under our assumptions. 

Note that the inequahty |x| ^ ka"^ is satisfied for d = 0{a) if fca ^ 1. 
Thus, formula (I22p is applicable in a wide region. 

Let us estimate Q asymptotically, as a — ?■ 0. 

Take cross product of (fTlll with A'' and integrate the resulting equation 
over S to get 

Io + h + l2 = 0, 
where Iq is defined in ([23|) (see below), Ii is defined in ([26]) . and 

/2 := - [ ^—[N, [iV, V X V X [ g(s,t)J(t)dt]]ds. 
Js *w/^o 

Let us estimate the order of Iq as a — > 0. One has 

lo = / ([N,Eo] - -^[N, [N,Vx Eo]])ds. (23) 

In what follows we keep only the main terms as o — )• 0, and denote by the 
sign ~ the terms equivalent up to the terms of higher order of smallness as 
a — )■ 0. One has 

/oo := / [N, Eo]ds = / V X Eodx ~ V x Eo{0)V, 
Js Jd 

where dx is the element of the volume, O € D is a point chosen as the origin, 
and V is the volume of D. Thus, /qo = 0{a'^). Denoting by \S\ the surface 
area of S, one obtains 

loi ■■=- f -T^[N, [N, V X Eo]]ds = ^rV x Eq, (24) 
Js «^/^0 ^W/io 

where 

T:=I-b, hraj -=7^ [ N^{t)Nj{t)dt, (25) 

Pi Js 
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b = {bmj) is a matrix which depends on the shape of S, and / := 6mj is 
the unit matrix. Since ( = 0{a~'^) one concludes that /qi = 0{a^~'^), so 
1-^00 1 1-^01 1 5 because /qo = 0{a^) as a — 0. Thus, 

h ^ hi = 0{a^-^), a ^ 0. 

Let us consider ly. 

h = J^ {J^[N, [Wgis,t), J{t)]]dt + ^)ds := hi + |. (26) 



One has 



in = J^ds J^dt(Vg{s,t)N{s) • J{t) - m^^)- 

It is well known that Jg ^^^^ds = — 5, where go := ^^|^_^ | . Since 

gis,t) -gois,t) = ^ + 0{\s -t\), \s - t\ ^0, 
47r 

and \s — t\ = 0{a), one concludes that if a ^ then 

/ Mfl^d. ~ -1 (27) 

and 

[ dt [ dsVg{s,t)N{s) ■ J{t) := /3Q, (28) 
Js Js 

where the matrix /3 is defined by the formula: 

P := iPmj) := / ^^Nj{s)ds. (29) 

Therefore, 

h^{I + P)Q. (30) 

Matrices (3 and h for spheres are calculated at the end of Section 1. 

Let us show that I2 = 0(a^~'^) and therefore I2 is negligible compared 
with /q as a — )• 0. If this is done, then equation for Q is 

Q = -^{I + a)TV xEq, (31) 

lljJjJLQ 

where 
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and the matrix I + /3 is invertible. From pip it follows that 

E{x) = Eo{x) - ^[V,g{x, O), {I + a)TV x Eq{0)\. (32) 

This equation is our first main result which gives an analytic formula for the 
solution of the EM wave scattering problem by a small body of an arbitrary 
shape, on the boundary of which an impedance boundary condition holds. 

In the far zone r := \x\ — > oo one has Vxg{x, O) = ikg{x, 0)x^ + 0(r~^), 
where x^ := x/r is a unit vector in the direction of x. Consequently, for 
r — > oo one can rewrite formula (j32p as 

E{x) = Eo{x) - C\S\(^) —[x'',{I + a)TVxEo{0)]. (33) 

This field is orthogonal to the radius-vector x in the far zone as it should 
be. 

Let us show that the term I2 is negligible as a ^ 0. Remember that 
curlcurl = graddiv — and 



-V^ / g{x,t)J{t)dt = / g{x,t)J{t)dt. 
Js Js 



Consequently, 



-iujfiol2^C / ds[N,[N, graddiv / g{x,t)J{t)dt]]\;j;^s ■= 
Js Js 

Since the function J(t) is assumed smooth, one has 

div [ g{x,t)J{t)dt = f g{x,t)divJ{t)dt, divJ 



'21- 



s Js 



dtr^ 



summation is understood here and below over the repeated indices, and 
divJ, where J is a tangential to S field, is the surface divergence. Further- 
more, 



grad I g{x,t)——dt = ep / g{x,t)———dt, 

Js *^^ni Js OtpOtm 

where the relation ^^g^'^^ = was used, and an integration by parts 

with respect to tp has ^een done over the closed surface S. Therefore 

h < c\C\\ [ ds [ \g{s,t)\ = 0{a^-^) « /q- 
Js Js 
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The constant c > here is a bound on the second derivatives of J on S. 
Example of calculation of matrices f3 and b. 

Let us calculate /3 and b for a sphere of radius a centered at the origin. 
It is quite easy to calculate b: 

If 1 

bjm = NjNrndt = -Sjm- 

\^\ JS "3 

Note that \S\ = Aira^, dt = sin 9d6d(j), and Nm is proportional to the 
spherical function Yi^m; so the above formula for j ^ m follows from the or- 
thogonality properties of the spherical functions, and for m = j this formula 
is a consequence of the normalization |A^| = 1. 
It is less simple to calculate matrix /?. One has 

47r7o Jo \s-t\-^ 

By symmetry one can choose t = (0,0,1) and let the z— axis pass through t 
from the origin. Then |t — = 2 — 2cos0, 1 — cos9 = 2sin^(^/2), and 

1 Z'^" r .„sin0fs™,5: 



Am = -r d(t> de 



^'^ 47ryo Jo 8sin3(0/2) 
This can be written as 

'"^ Ida h Jo siii'(9/2) 

If j 7^ m, then f3jm = 0, as one can check. If j = m < 2, then 5^3 = and 
f3jj = ^, as one can check. Finally, 

f27r 



1 f r ,m /n^ COS f - COS^ _ 1 

sin2(0/2) ~ 6' 



hi = ^ r d<^ rde cos{e/2) 



Thus, 



1 

6' 

Let us formulate the result of this Section as a theorem. 

Theorem 1. If ka <^ 1, then the solution to the scattering problem 
dH-dni) is given by formula (j32]) . 

In the next Section the EM wave scattering problem is studied in the 
case of many small bodies (particles) whose physical properties are described 
by the boundary impedance conditions. 
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3 Many-body EM wave scattering problem 



Consider now EM wave scattering by many small bodies (particles) Drn^ 
1 < m < M = M{a), distributed in an arbitrary bounded domain 17. Here 
a is the characteristic dimension of a particle, M{a) is specified below, and 
Xm G Dm is a point. Assume for simplicity of formulation of the problem 
that the particles are of the same shape and orientation, and their physical 
properties are described by the boundary impedance p6p . where h = h[xm), 
and h[x), Re/i(x) > 0, is a continuous function in an arbitrary large but 
finite domain 17 where the particles are distributed. For simplicity let us 
assume that the particles are distributed in a homogeneous medium with 
parameters eo,//o- 

The distribution of small particles is given by the following law: 

A/'(A) = a-^^-'') [ N{x)dx{l + o{l)), (34) 

where A is an arbitrary open subset of 17, N{x) > is a given function 
which one can choose as one wishes, A/'(A) is the number of particles in 
A, and n S (0,1) is a parameter from formula (I16p . The distance d = 
d{a) ^ a between the neighboring particles is assumed to satisfy the relation 
\mia-^Q ad"^ {a) = 0. The order of magnitude of d as a — )■ can be found 
by a simple argument: if one assumes that the particles are packed so that 
between each neighboring pair of particles the distance is d, then in a unit 
cube Ai the number of particles is 0{d~^), and this amount should be equal 
to A/'(Ai) = 0(a-(2-«)). Therefore 

d = 0(a(2-«)/3). 

Denote D := U^=i D^, D' := \ D. 

The scattering problem consists of finding the solution to 

V X V X = k'^E in D' , (35) 

[N, [E, N]] - Cm[N, VxE] on Sm, 1 < m < M, (36) 

E = Eo + VE, (37) 

where Eq is the incident field, which satisfies equation (I35p in M^, and ve 
satisfies the radiation condition. Existence and uniqueness of the solution 
to the scattering problem (f35]l - (f37|l is known. Let us look for the solution of 
the form 

M 

E{x) = Eo{x) + Y,'^ ^ / 9{x,t)Jm{t)dt. (38) 

m=l 
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Define the effective (self-consistent) field, acting on the j—th particle, by the 
formula 

M 

Ee{x) = Eoix)+ g{x,t)Jm{t)dt. (39) 

Because of our assumption d ^ a, any single particle Dj can be considered 
as being placed in the incident field Ee{x). Therefore 

*^ C 151 

Ee{x) = Eo{x)- V ^^[V xgi^X, Xfn ),(/ + a)rV X Se(x„^)]. (40) 

Note that C,\S\ = h{xm)csa?~'^ , where the constant cs depends on the shape 
of S and is defined by the formula 

Thus, formula ()iO]) implies 

M 

Ee{x) = Eo{x) - ^^a"^-"" V h{x^)[Vxg{x,Xm),{I + a)TVxE^{x„,)]. 

(41) 

Denote Am := V x Ee{xm)- Applying operator Vx to ([IT]) and then setting 
X = Xj, one obtains a linear algebraic system (LAS) for the unknown A^,'- 

M 

Aj = Aoj - a^"'' ^ h{xmWx X [V xg{x,Xm)\Al + Ol)'rAm]\x=Xj, 

lUJflQ . 

(42) 

where 1 < j < M. If A^ are found, then formula (j40p yields the solution to 
many-body EM wave scattering problem. The magnetic field is calculated 
by the formula H{x) = -^^j^^, if E{x) is found. 

Let us now pass to the limit in equation (|41|) as a 0. 

Consider a partition of into a union of cubes Ap, 1 < p < P. These 
cubes have no common interior points, and their side i ^ d, I = l{a) — > 
as a — )■ 0. Choose a point Xp € Ap, for example, let Xp be the center of Ap. 
Rewrite ()4ip as 

p 

£;e(xg) = £^o(a;<7) - /i(xp)[Va;5(x,Xp)|a;=j: , (/ + a)rV x ^e(a;p)] 



a' 



(43) 
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It follows from equation that 

-.2-K 1 ~ iV(xp)|Ap|, a^O, 



a 



where |Ap| denotes the volume of the cube Ap. Thus, ()43p can be written 
as: 

y~] h{xp)N{xp)[S/xgix,Xp)\x=x„,{I + Oi)TV X Ee{xp)]\Ap\. ^ 

This equation is a discretized version of the integral equation: 

E{x) = Eoix) - X / g{x, y){I + a)rV x E{y)N{y)h{y)dy. (45) 

Equation (j45|) is the equation for the limiting, as a — )• 0, effective field in the 
medium in which M(o) small impedance particles are embedded according 
to the distribution law 

Let us summarize our result. 

Theorem 2. If small particles are distributed according to ()34p and 
their boundary impedances are defined in (jl6p . i/ien the solution to a many- 
body EM wave scattering problem (|35]) - (i37|l is given by formula (jiOl) . where 
the quantities V x E{xm) '■= o-Te found from the linear algebraic system 
(jl2|) . and i/ie limiting, as a ^ 0, electric field E{x) in solves integral 
equation ([I5]) . 

Let us discuss the novel physical consequences of our theory. Apply- 
ing the operator curlcurl to equation (I45p and taking into account that 
curlcurl = graddiv — A, divcurl = 0, and —Ag{x, y) = k'^g{x, y) + 5{x — y) 
one obtains 

V xV X E = k^E - V X (-^U + a)rV x E(x)N(x)h(x)) . (46) 

\iujfj,o ' 

The novel term, which is due to the limiting distribution of the small parti- 
cles, is 

T := -V X (-^{1 + a)rV x E(x)N{x)h{x) 
\iujfio 

Let us make a simple assumption that N{x)h{x) = const. Physically this 
means that the distribution of the small particles is uniform in the space 
and the boundary impedances of the particles do not vary in the space. 
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Denote -j^^Nh := ci. Assume also that a and r are proportional to the 
identity matrix. This, for example, happens if all Sm are spheres. Then 
T = — C2V X V X E{x), where C2 = const. 

In this case the novel term T can be interpreted physically in a simple 
way: it results to a change in the refraction coefficient of the medium in il. 

Indeed, by taking the term T to the left side of ()46p and dividing both 
sides of this equation by 1 + C2 , one sees that the coefficient k'^ is replaced by 
^1 •= 2+C2 • Si'^*^^ is a complex number, the new medium is absorptive. 
The assumption Re^ < implies that Im/c^ > 0. Since h{x) and N{x) can 
be chosen by the experimentalist as he wishes, their dependence on x can be 
chosen as he wishes. Therefore, distributing many small particles in one 
can change the refraction coefficient of the medium in a desirable direction. 

4 Conclusions 

The main results of this paper are: 

1. Equations ()32p and ()31|) for the solution of the EM wave scattering 
problem ([J])-® by one small body of an arbitrary shape. 

2. A numerical method, based on linear algebraic system (LAS) ()42p . see 
formulas (|4ip and (j42p . for solving many-body EM wave scattering problem 
in the case of small bodies of an arbitrary shape. 

3. Equation (j45p for the electric field in the limiting medium obtained 
by embedding M(a) small impedance particles of an arbitrary shape, dis- 
tributed according to ([M|) . as a — > 0. 
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